THE  SINGULAR  SURFACE  IN  THE  PARAMETER 
SPACE;  APPLICATIONS 


Paul  Franklin  Van  Tassel 


DUDLEY  KNOX  LIBRARY  ^, 

NAVAL  WSTGRAJU  ATE  ^HOOI, 

MONTEREY,  CALIFORNIA  93SW 


AL  POSTGRADUATE  SCH 

Monterey,  California 


L 


*sam*?    fit    ^*»^ 


THE 

SINGULAR 

SURFACE  IN  THE  PARAMETER 

SPACE ;   APPLICATIONS 

by 

Paul 

Franklin  Van  Tassel 
December  1975 

Thesis 

.  Advisor 

:                  G.J. 

Thaler 

Approved  for  public  release;  distribution  unlimited. 


T17217Z 


SECURITY   CLASSIFICATION   OF    THIS  PAGE  (When  Data  Entered) 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER 


2.  GOVT  ACCESSION  NO 


4.      TITLE   (end  Subtitle) 


The  Singular  Surface  in  the 
Parameter  space;  Applications 


7.  AUTHORS 


Paul  Franklin  Van  Tassel 


READ  INSTRUCTIONS 
BEFORE  COMPLETING  FORM 


3.      RECIPIENT'S  CATALOG   NUMBER 


5.     TYPE  OF  REPORT  a  PERIOD  COVERED 

Master's  Thesis 
December  1975 


6.  PERFORMING  ORG.  REPORT  NUMBER 


8.  CONTRACT  OR  GRANT  NLM8ER(.; 


9.  PERFORMING  ORGANIZATION  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93943 


10.  PROGRAM  ELEMENT.  PFOJECT  TASK 
AREA  a  WORK  UNIT  NUMBERS  ' 


II.  CONTROLLING  OFFICE  NAME  AND  ADDRESS 

Naval  Postgraduate  School 
Monterey,  California  93940 


12.  REPORT  DATE 

December  1975 


13.  NUMBER  OF  PAGES 

87 


U.     MONITORING  AGENCY  NAME  a    AODRESSfi/  different  from  Controlling  Otllce) 


Naval  Postgraduate  School 
Monterey,  California  9  3940 


IS.     SECURITY  CLASS,  (ol  Mia  riport) 


Unclassified 


15a.     DECL  ASSIFI  CATION/  DOWN  GRADING 
SCHEDULE 


16.     DISTRIBUTION   ST  ATEMEN  T  (of  thlt  Report) 


Approved  for  public  release;  distibution  unlimited 


17.     DISTRIBUTION  STATEMENT  (of  the  abetract  entered  In  Block  30,  It  different  from  Report) 


18.     SUPPLEMENTARY  NOTES 


19.     KEY  WORDS  (Continue  on  reveree  elde  It  neceeeary  end  Identity  by  block  number) 

Parameter  Space,   Parameter  Plane,   Singular  Surface, 
Singular  Line 


20.     ABSTRACT  (Continue  on  reveree  elde  It  neceeemry  mid  Identity  by  block  number) 


The  basic  theory  of  the  parameter  plane  and  the 
singular  line  is  revalued.   Previously  unknown  properties  of 
the  singular  line  and  its  relationship  to  the  parameter 
plane  are  introduced. 


DD  ,£"»  1473 

(Page  1) 


EDITION  OF    1  NOV  85  IS  OBSOLETE 

S/N   0102-014-6601   | 


SECURITY  CLASSIFICATION  OF  THIS  PAGE  (When  Dele  Kntered) 


CtiCU'flTY   CLASSIFICATION   OF    THIS  PAGECWTian  Del«  Entmfdj 

The  singular  surface,  a  hypersurface  of  singular  value 
in  the  parameter  plane,  is  introduced. 

Compensation  techniques  using  the  singular  line  are 
demonstrated  and  a  computer  program  for  determining  the 
existance  of  singular  lines  is  presented. 


DD  Form   1473     (BACK) 

.  1  Jan  73  , M 

S/N      0102-014-6601  SECURITY  CLASSIFICATION  OF  THIS  PAGEfWh.n  D.f.  Ent.fd) 


The  Singular  Surface  in  the  Parameter  Space;  Applications 


by 


Paul  Franklin  Van  Tassel 

-   4 

Lieutenant,  United  States  Navy 

BSEE,  Purdue  Univeristy,  1S70 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  ELECTRICAL  ENGINEERING 

from  the 


NAVAL  POSTGRADUATE  SCHOOL 
December  1975 


-V  KNOX  LIBRARY 

L  POSTGRADUATE  SCHOOL 

-REY.  CALIFORNIA  93S40 


ABSTRACT 

The  basic  theory  of  the  parameter  plane  and  the 
singular  line  is  reviewed.  Previously  unknown  properties  of 
the  singular  line  and  its  relationship  to  the  parameter 
plane  are  introduced. 

The  singular  surface,  a  hypersurface  of  singular  values 
in  the  parameter  plane,  is  introduced. 

Compensation  techniques  using  the  singular  line  are 
demonstated  and  a  computer  program  for  determining  the 
existance  of  singular  lines  is  presented. 


TABLE  OF  CONTENTS 

I.  INTRODUCTION 10 

II.  THE  SINGULAR  LINE  IK  THE  PARAMETER  PLANE 13 

A.  THE  PARAMETER  PLANE 13 

B.  THE  SINGULAR  LINE 23 

III.  PROPERTIES  OF  THE  SINGULAR  LINE 29 

IV.  THE  SINGULAR  SURFACE  IN  PARAMETER  SPACE 47 

A.  PARAMETER  SPACE 47 

B.  SINGULAR  SURFACES 50 

V.  COMPENSATION  TECHNIQUES  USING  THE  SINGULAR  LINE 54 

VI.  CONCIUSICNS  AND  RECOMMENDATIONS  FOR  FURTHER  STUDY-  61 

APPENDIX  A  PARAM  A ■—  63 

APPENDIX  E  SINGULAR  LINE 64 

COMPUTER  PROGRAM  PARAM  A 66 

COMPUTER  PROGRAM  SINGULAR  LINE 74 

BIBLIOGRAPHY 86 

INITIAL  DISTRIBUTION  LIST 37 


LIST  OF  TABLES 


TABLE 


PAGE 


II- 1 


POINTS  FROM  FIGURE  2.2 21 


II-2    POINTS  FROM  FIGURE  2.3 28 


II-3    POINTS  FROM  FIGURE  2.3 28 


LIST  OF  FIGURES 

FIGURE  PAGE 

2.1  A  POINT  IN  THE  S-PLANE 14 

2.2  PARAMETER  PLANE  EXAMPLE,  FOURTH  ORDER  POLYNOMIAL. 22 

2.3  SINGULAR  LINES  ON  THE  PARAMETER  PLANE, 

FOURTH  ORDER  POLYNOMIAL 27 

3.1  AIRCRAFT  PITCH  STABILIZATION  SYSTEM 30 

3.2  AIRCRAFT  PITCH  STABILIZATION  PARAMETER 

PLANE  WITH  ZETA=0. 6-OMSGA= 12. 74  SINGULAR  LINE. ..33 

3.3  PARAMETER  PLANE  FOR  PITCH  STABILIZATION 

SYSTEM,  ZETA=0.6 37 

3.4  SINGULAR  LINE  THROUGH  THE  PARAMETER  PLANE, 

FOURTH  ORDER  POLYNOMIAL 40 

3.5  SINGULAR  LINE  IN  THE  PARAMETER  PLANE,  EXAMPLE  V..43 

3.6  SINGULAR  LINES  IN  THE  PARAMETER  PLANE, 

EXAMPLE  V  -  EXPANDED 44 

4.1  THIRD  ORDER  SYSTEM  WITH  DUAL  CASCADE 

COMPENSATION 4  7 

4.2  PARAMETER  SPACE  WITH  ZETA=0.5  AND  OMEGA=1.7 

SURFACES 4  9 

4.3  SINGULAR  SURFACE  MODEL 51 

4.4  SINGULAR  SURFACE  DRAWING 52 

5.1    AIRCRAFT  PITCH  STABILIZATION  SYSTEM 55 


5.2  SINGULAR  LINES  FOR  AIRCRAFT  PITCH  STABILIZATION 

SYSTEM 5  8 

5.3  STEP  RESPONSE  FOR  POINT  A 1.59 

5.4  STEP  RESPONSE  FOR  POINT  B 59 

5.5  AIRCRAFT  PITCH  STABILIZATION  SYSTEM 

PARAMETER  PLANE 60 


ACKNOWLEDGEMENT 


The  author  wishes  to  express  his  sincere  appreciation 
to  Dr.  George  J.  Thaler,  Professor  of  Electrical 
Engineering,  Naval  Postgraduate  School  for  his  aid  in  the 
selection  of  this  subject  and  his  assistance,  encouragement 
and  guidance  in  the  preparation  of  this  work.  Appreciation 
is  also  extented  to  the  W.R.  Church  Computer  Center  for  the 
cooperation  and  assistance  rendered  in  the  preparation  of 
this  document. 


I.   INTRODUCTION 


Compensation  of  linear  control  systems  by  the  classical 
methods  cf  root  locus  and  frequpncy  response  is  adequate 
when  there  is  a  single  variable.  However  when  there  are  two 
or  more  variables  compensation  of  physical  systems  for 
desired  system  dynamics  becomes  a  complex  task  [1]. 

In  modern  control  engineering  systems  characteristic 
equations  with  two  or  more  unknown  parameters  result.  In 
1959,  Kitrovic's  Method  [2]  was  developed  to  portray  the 
relationship  between  the  roots  of  the  characteric  equation 
and  two  coefficients  of  the  system.  This  method  allowed  the 
two  lowest  ordered  coefficients  of  the  characteristic 
equation  to  be  varied.  Siljak  expanded  this  method  to  tae 
Coefficient  Plane  to  allow  the  variation  of  any  two  of  the 
coefficients  of  the  characteristic  equation.   For  example: 


n 
i      j    r-i    k 
f(s)=Bs   +5s  +  )   as  (1-1) 

i      j      t—»      k 

k=0 
k*i,j 

where:     a   are  real  constants 
k 

B   and  B   are  real  variables 


Application  of  this  method  resulted  in  constant  zeta,  omega, 
and   sigma  curves  on  the  B  -B   plane  specifying  the  roots  of 

equation  (1-1)  for  any  choice  of  £  and  CO  • 

n 


In   1964,   Siljak   [3]  extended  Mitrovic's  work  so  that 
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the  two  variable  parameters  may  appear  linearly  in  the 
coefficients  of  the  characteristic  equation.  consider  the 
linear  case  of  the  characteristic  equation: 


n 
f(s)  =  J]    a^s   =  0  (T-2) 

k=0 

where:  a   =  bCl  +  c  $  +  d  (1-3) 

k    k     k     k 

and:  and    are  real  variable  parameters 


Application  of  this  method  resulted  in  constant  zeta,  omega 
and  sigma  curves  for  any  choice  of  Q,  and  ft  on  the  &- Q  plane 
specifying  the  roots  of  equation  (1-2) . 

Siljak*s  Parameter  Plane  Method  provides  the  designer 
with  a  simple  procedure  for  factoring  a  characteristic 
equation  and  displaying  the  results  in  a  parameter  plane 
diagram.  It  allows  the  designer  to  obtain  information  about 
system  stability  and  the  affect  of  parameter  variations  and 
adjustments  on  stability.  In  the  general  case  the  Parameter 
Plane  Method  can  be  applied  to  any  control  problem  in  which 
it  is  necessary  to  determine  how  the  variation  of  one  or 
more  parameters  effect  the  root  location  of  a  system. 

In  1967,  Bowie  [4]  discovered  that  a  limitation  existed 

in   parameter  plane   theory.    While   attempting  to  solve  a 

sixth  order  characteristic  equation,   Bowie   found   that   he 

could   not   find   a   complete  set  of  roots  in  an  area  of  the 

parameter  plane  that  should  have  contained  only  real   roots. 

By    substituting    the    values    of   Q,   and  Q        into    thu 

characteristic   equation   he   found   the   roots.     It    was 

discovered  that  a  set  of  complex  roots  existed  corresponding 

to   a  C -  CO   pair   in   the   left-half    s-plane.     Further 
n 
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investigation  showed   that  this  Z,~U)    pair  formed  a  straight 

n 

line  in  the  parameter  plane.   This  was  contrary  to  the  then 

existing   parameter   plane   theory   which  stated  that  a  £-CU 

n 

pair  will  only  occupy  a  point  in  the  parameter  plane. 


Bowie   called   this   line   of  constant  £  -constant^  the 

n 

SINGULAR  LINE  as  it  was  formed  by  a  singular  value  of 

C      and  CO   The   singular   line  added  a  new  dimension  to  the 
n 

parameter  plane  because  it  allows  the  designer  to   hold   one 

root   location   constant   while   moving   others   to   achieve 

desired  system  characteristics  by  varying  (L  and  0  . 
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XI.   THE  SINGULAR  LINE  IN  THE  PARAMETER  PLANE 

A.   THE  PARAMETER  PLANE 

The  Parameter  Plane  Method  was  developed  by  Siljak  [3] 
as  an  extension  to  Mitrovic's  Coefficient  Plane.  The 
parameter  plane  is  a  mapping  of  points  in  the  s-plane  into 
the  G<  -  G>  plane.  The  following  is  the  mathematical 
development  of  that  transformation. 

Consider  the  characteristic  equation: 


f  (s) 


11 

■r 

k=0 


(2-1) 


where : 


and: 


a   (k  =  G, 1 , 2 , . . . , n)  are  real  variables 
k 


S  =  -ct+Jw  =  - £w  ±jw   \/   1  — I 

n   n 


(2-2) 


Examination  of  figure  2.1  shows 


-  £  =  Cos# 


and 


V 


By   substituting   equation 
expressed  as: 


-  'C      =  Sinfl 
(2-3) ,   equation   (2-2) 


(2-3) 


mav      oe 


and 


s   =  CO      {Cos$      ♦    jSinfl) 
n 


s      =  CO      (Cosk0   +    jsinkfl) 

D 


(2-4) 
(2-5) 
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♦  Im 


-P.e  4- 


Figurc  2.1 
A  POINT  IN  THE  S-PLANE 
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Since 


-1 
0  =  Cos   (-£)  (2-6) 


it  follows  that 

Coskfl  =  Cosk[Cos(-0]  (2-7) 

and  Sink0  =  Sink[  Cos  (-£)  ]  (2-8) 

Equation  (2-7)  can  be  expressed  as  the  Chebyshev  function  of 

the   first   kind,  I  (-£)  and  similarly  equation  (2-8)  can  be 

k 

expressed  in  terms  of  the  Chebshev  function   of   the   second 

kind,   Sin   U  £-  \   ].   Substituting  equations  (2-7)  and  (2-8) 
k 

equation  (2-5)  becomes 


k     k 
s   =  a/  [T  (-/-)  +  jSin  U  (-£)  ]  (2-9) 

n   k  k 


k 
where:     T  (-£)  f=  (-1J   I  {£) 


k+1 
and  U  <-£>  =  (-1)     (J  (C) 

k  k 


(2-10) 


Substituting  equation  (2-10)  into  equation  (2-9),  then 
substituting  equation  (2-9)  into  equation  (2-1),  equation 
(2-1)  becomes 


n  n 

V1    k    k  c— i    k    k+1 
f(s)=)  a  W(-1)  T  ft-)  +  jsin  >  a  W  (-1)    U  (J)  =0  (2-11) 

Z_j  k  n      k  i— t    k  n        k 

k=0  k=0 

For  any  non-zero  value  of  both   the   real   and   imaginary 

parts   of   equation   (2-11)  must   go   to  zero  independently 
therefore, 
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n 


Ek    k 
a,   (-1)  T, 


k  =  0 


k  n 


(2x12) 


u 

I 

k=0 


k+1 

k  k        k 


The   Chebyshev   functions   of  the  first  and  second  kind 
are  given  by  the  recursion  formulae 


T    (*)  =  2T  (J)  -  T    (Y) 

k+1        k       k-1 


0    ft)  =  2U<  )  -  0    (J) 

k+1  k-1 


(2-13) 


where:     T  (t)  =  1 
0 


u    or.)  =  -1 

-1 

Oft)  =  o 

0 

a  (5)  =  1 


To  eliminate  T  (})  from  equation  (2-12)  it  is  noted  that 
k 


VT)  =   \™  'Vi(^ 


(2-14) 


then  equation  (2-12)  may  be  rewritten 


k        k    k-1 


t-"1  k        Y.  r~»  k        k    k-1 

)     (-1)    aCOC    U     <T)]->        (-1)    aWU         (J)    =0  (2-15a) 

£__>  k   n        k  i— j  k   n 


k=0 


k=0 


k+1        k 
(-1)  a   COO    O) 

k   n   k 


=   0 


(2-15b) 


Examination  of  equation  (2- 15a)  shows  that  its  first  part  is 
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a  linearly  dependent  solution  to  equation   (2-15b)   for   ail 
D^n-zero  values  of  zeta  therefore  equation  (2-15)  becomes 


t 

k=0 


k   k 
k  n  k-1 


(2-16) 


Ek   k 
(-1)  a    UK)     =  0 
k  n  k 


k=0 


For  the  linear  case: 


a   =  CLb   +  Sc   +  d 
k     k   r  k    k 


(2-17) 


where: 
and : 


b  ,  c  ,  and  d   are  real  constants 
k   k       k 

CL-and  ^  are  the  variable  parameter: 


Substitution  of  equation  (2-17)  into  equation  (2-16)  yields: 


B<X    +    C  Q>   *    D 
1  V  1 


=    0 


(2-13) 


and 


B  <X  +    C    fl    +    D 
2  2V  2 


=    0 


where: 


\ — '  k        k 

1  =2_j~1)  b-  -U] 


k        k 


k  n   k-1 


(-D    b      urn 

k   a   k   J 


k=0 


k=0 


Ek        k 


k  =  0 


k   n   k-1 


c.    =Y    (-1)    c        &,  (I)        <2~19) 
2      /_»  k   n    k 
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n 

k   k 


k  n  k-  I 


k=0 


Application  of  Cramer's  Rule  for  the  solution  of 
simultanous  linear  equations  to  equation  (2-18)  yields  the 
desired  parameter  plane  solution  equation: 


CD-CD  BD-3D 

12     2  1  ^21     12 

£  = S= (2-20) 

BC-BC  *    BC-BC 
12     2  1  12     2  1 


Siljak  [3]  defined  the  parameter  plane  as  a  rectangular 

plot  with    as  the  abscissa  and    as  the  ordinate.   Equation 

(2-20)  gives  Q,  and  ft  as  a  function  of  CO  and  \   .   By  f ixina  UJ=a* 

>  n  i 

and  varying  (0  over  the   range   from   zero   to   infinity   in 
n 

equation   (2-20),   a   curve  Is  produced  in  the  cC-@  plane  for 

the  fixed  value  of  zeta.   This   curve   specifies   the   «^  and^> 

pairs   which   will   cause   equation   (2-1)  to  have  a  pair  of 

roots  corresponding  to  the  fixed  value  of   zeta.    Similarly 

by   fixing  CO   =  CO    an(*   varying   zeta   from  -1  to  +1  in 
n       ni 

equation  (2-20),  a  curve  is  produced  in  the  CL - Q    plane   for 

the   fixed   value   of  CO   •    This  curve  specifies  the  Q,   and  $ 

ni 

pairs  which  cause  equation  (2-1)  to  have   a   pair   of   roots 

correspondmq   to   the   fixed   value   of  ^j  .   Thus  equation 

n 

(2-20)  allows  the  mapping  of  complex  points  in   the   s-plane 

to  the  parameter  plane. 

To  transform  the  real  axis  points  in  the  s-plane  into 
the  parameter  plane,  replaces  s  in  equation  (2-1)  by  -«■  to 
obtain: 


18 


n 

f(s>  =  2^  v~CT)   =  °  (2~21> 

k-0 

Now  equation  (2-17)  is  substituted  into  equation   (2-21)   to 
obtain: 


n 

k 


Y  ( bj^+^ck  +  d>  (-a)  =  °  <2"22) 


k=0 
Simplifying  equation  (2-22)  obtains: 


B(<0   +^C(a)  +  D(<r)  =0  (2-23) 


where 


n 

k    k 
b  ■  " 

k=0 


n 

k    k 

C(a)      =       )  C  a 


ii 

dm  =  r 


n 

k   k 

(-1)   d  cr 

k 
k=0 


For   a   given  value  of  <r  the  functions  B(<0  ,  C(cr)  ,  and  D(o-) 
are  constants  and  equation  (2-20)  results  in  a  straight  line 
in   the   parameter  plane  with  a  locus  of  point  corresponding 
to  the  real  roots  s=-a. 

A  digital  computer  program  PARAM  A  was  written  by 
Nutting  [5]  to  aid  in  solving  the  characteristic  equation 
and  plotting  the  constant  zeta,  omega,  found  at  the  end  of 
this  thesis  and  a  brief  explanation  is  found  in  Appendix  A. 
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EXAMPLE  I: 

Consider  a  system  whose  characteristic  equation  is   the 
fourth  order  polynomial: 

4      3         2 
f(s)=s+8.5s  +(5a+28)s  +  ( 1  2  .  5GL+25f  +  42  .  5)  s  +  50£  +  25  =  0     2-25) 

Figure   2.2   shows   the   parameter   plane   for  the  following 
values  of  zeta,  omega,  and  sigma: 


OMEGA 


ZETA 


SIGMA 


3.0 

0.00 

0.5 

4.0 

0.30 

1.0 

5,0 

0.35     . 

2.0 

6.0 

0.40 

3.0 

7.0 

0.45 

4.0 

8.0 

0.50 

5.0 

9.0 

1.00 

6.0 

Examination   of   figure   2.2   shows   the  region  of  stability 
bounded  by  ^=0.0,  *J  =  1.0,  and(T~0.0.   All  values  ofcC  and  ^ 

inside  rhis  region  will  result  in  a  stable  system.  The 
roots  of  the  system  can  be  determined  by  selecting  a  point 
in  the  plane  and  reading  the  corresponding  values  of  zeta, 
omega,  and  sigma.  Zeta  and  omega  will  form  a  complex  of 
roots  so  that: 


s  =  -JCu  ±  jcj 
n      n 


rVT 


and  the  sigma  poii>t  will  form  the  real  roots  of  the  form: 


S     =     -<r 


Although  it  may  not  be  immediatly  obvious,  there  are 
two  sigma  points  through  each  point  in  the  region.  Careful 
interpolation   or  the  plotting  of  additionalines  is  reguired 
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to  obtain  the  exact  roots  of  the  system.  But  this  is 
generally  unnecessary  as  the  designer  will  interested  in  the 
general  location  of  a  specific  pair  of  roots. 

Four  points  have  been  selected  from  figure  2.2.  the 
position  of  these  points  was  measured  with  an  engineer's 
scale  and  the  resulting  value  of  and  were  substituted  in 
equation  (2-25).  The  roots  of  the  equation  were  then  solved 
by  the  use  of  a  polynomial  rootfinder  program  on  tne  IBM 
360/67.  Table  II- 1  shows  the  coordinates  of  the  points,  the 
resulting  roots  and  the  corresponding  values  of  zeta,  omega, 
and  sigma.  As  can  be  seen  the  parameter  plane  provides  a 
simple  means  of  determining  a  root  location  while  variable 
parameters  are  adjusted. 


17>  71   T3  T   T?        T  T—  I 
1  nt>*ju      -A-  JL       I 


POINT 

ALPHA 

BETA 

COMPLEX  HOOTS 

REAL  ROOTS 

OMEGA 

ZETA 

A 

6.  1 

4.0 

-l.461±j5.814 

-1.557 
-4.020 

5.9  95 

.244 

B 

6.  1 

3.0 

-1.815±j5.78  3 

-1.076 
-3.794 

6  .  0  6  1 

.299 

C 

7.0 

1.975 

-2.278±j6.075 

-0.998 
-2.945 

6.4  90 

.351 

D 

2.8 

0.725 

-2.308±j3.982 

-1.004 
-2.881 

4.603 

.501 

POINTS    FROM    FIGURE    2.2 
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a-scale:  4  units/inch 
/?-scale:  1  unit/inch 


Figure  2.2 


PARAMETER  PLANE  EXAMPLE 
FOURTH  ORDER  POLYNOMIAL 
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3.   THE  SINGULA*  LINE 

As  stated  in  Chapter  I  Bowie  [4]  discovered  the 
existance  of  the  singular  line  in  the  parameter  plane.  In 
investagating  the  cause  of  singular  lines  Bowie  found  that 
all  previous  work  in  the  parameter  plane  was  correct  for  as 
far  as  it  went.  It  did  not,  however,  give  a  full 
explanation  of  the  case  when  both  the  numerator  and 
denominator  of  equation  (2-20)  formed  singular  matrices. 
Siljak  [3]  did  address  the  fact  that  this  condition  could 
exist,  nut  the  idea  of  a  singular  line  in  the  parameter 
plane  had  not,  up  to  that  time,  been  conceived. 

Bowie  found  that  in  order  for  singular  lines  to  exist 
in  the  parameter  plane  the  following  conditions  must  b<=  met: 


3  C   -  B    C   =0 
12     2  1 


CD   -CD   =0  (2-26) 

12     2  1 


B  D   -  B  D   =  0 
2  1     12 


where 


r— i     k   k  r— »     k   k 

1  ii-i.  k  ns  k- 1  *s  2      L-t  k  ns  k  -*s 

k=0  k=0 

n  i  n 

s  =E  <-A  lv,(v  c2  -E«-i»\  Lw  (2-27) 


k=0  k=0 

n  n 


Ek   k      w  r-i     k    k    v 

(-1)  a   u   ft  )     d  =)  (-1)  a   u  CI  ) 
k  ns  k-1   s       2  Is  k  ns  k   s 


.k=0  k=0 
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ind:     CO       ana"  X      are  tiie   constantto-constantt      values. 
ns  Js  s 


equations       (^-18)       are      solyed    for         in   teras    of  Q  the 
It    is    two   linearly    independent   equations: 


D  B  D  B 

n  1  1  2  2 

~  c  c  c         c 

1         1 


(2-28) 


Equation  (2-28)  is   the   singular   line   solution.  It 

s  a  straight  line  in  the  parameter  plane  with  a  slope  of 

C   and  it  intercepts  the   -axis  at  -D  /C  .  All  points  on 
1   "                                11 

line  represent  a  singular  value  of  zeta  and  a  singular 

e  of  omega  and  therefore   map  into   a   single   pair  of 
lex  points  in  the  s-plane  where: 


s  =  -X  CO   ±  jco 

s  ns      ns 


i-r 


(2-29) 


It   is  obvious  that  singular  lines  do  not  exist  for  all 

acteristic  equations  and  when  they  do   exist   they   will 

;t   only  for  the    -CO   pairs  which  will  satisfy  eguation 

s   ns 

j>6)  .   To  aid  in  determining   the   existance   of   singular 

is      and  plotting  them  in  the  parameter  plane  the  computer 

:jram  SINGULAR  LINE  has  been  written. 


Briefly,   the   program   uses   the  real  constants  of  the 

fficients  of  the  characteristic  equation  and  the   desired 

■ie      of  V   to  solve  for  value  for  any  real  positive  values 
s 

0   in  equation  (2-26)  .   The  Y   and  CO    values   are   then 
ns  s      ns 

<d  to  generate  equations  (2-27)  and  these  results  are  used 

test  for  the  existance  of  singular   lines.    If   singular 
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region  there  will  be  two  real  roots  and  a  pair  of  complex 
roots  corresponding  to  the  constant  zeta-constant  omega 
values  of  the  singular  lines. 

A   pair   of   points   nas  been  selected  along  two  of  the 

singular  lines.   Points  A  and  B  are  on  the  line  for  ^J  =0.7, 

CO  =13.28  and  point  C  and  D  are  on  the  line  for  3T  =  0 . 4 
n 

60=12.309.   Table  II-2  shows  the  coordinates,  the  resulting 
n 

roots   and  the  values  of  zeta  and  omega.   For  Table  II-2  the 

coordinates  of  the  point  were  measured   with   an   engineer's 

scale  and  the  values  of  alpha  and  beta  were  substituted  into 

eguation   (2-30)   and   the   roots   were   computed   using    a 

polynomial  rcotfinder  subroutine  on  the  IBM  360/67. 

Table  II-3  is  for  the  same  points,  tne  only  difference 
being  that  the  coordinates  were  computed  using  the 
slope-intercept  method  utilizing  the  printed  output  from  the 
program  SINGULAR  LINE. 

Table  II-3  shows  that  this  method  of  obtaining  system 
roots  provides  an  accurate  method  for  obtaining  singular 
lines  and  that  the  lines  a  for  constant  zeta-constant  omega 
values,  while  Table  II-2  shows  that  reasonably  accurate 
results  can  be  obtained  using  a  minimum  of  equipment. 


26 


«-scale:  5  units/inch 
i8-scale:  1  unit/inch 


Figure  2.3 


SINGULAR  LINES  ON  THE  PARAMETER  PLANE 
FOURTH  ORDER  POLYNOMIAL 
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TABLE  II-2 


POINT 

ALPHA 

BETA 

COMPLEX    ROOTS 

REAL    ROOTS 

OMEGA 

ZETA 

A 
B 
C 
D 

20. 
10. 
10. 
15. 

19.67 

10.  5 

5.667 

8.  167 

-9.244±j10.87 
-9.281±j9.670 
-4.839±j11. 19 
-4.91S±j11.27 

-.9326 

-422.6 

-1.051 
-222.4 

-.9971 

-231.3 

-1.045 
-326 . 7 

13.40 
13.40 
12.  20 

12.30 

.639 
.693 
.39"? 

.4  00 

P 

OI1JTS    FROM    FIGUI 

IE    2.3 

TABLE  II-3 


point 

ALPHA 

BETA 

COMPLEX    ROOTS 

REAL    ROOTS 

OMEGA 

ZETA 

A 

20. 

t9.S6 

-9.298±j9.466 

-1.071 

-422. 3 

13.40 

.700 

B 

10. 

10.52 

-9.298±j9.486 

-1.051 
-222.4 

13.40 

.700 

c 

10. 

5.755 

-4.920±j11.27 

-.9982 
-  23  1 .  2 

12.30 

.400 

D 

15. 

8.225 

-4.920±j11.27 

-1.011 

-331.2 

i — , 1 

12.30 

.4  00 

POINTS  FROM  FIGURE  2.3 
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III.   PROPERTIES  OF  THE  SINGULAR  LINE 


To  fully  understand  the  relationship  of  singular  lines 
to  the  parameter  plane  it  is  necessary  to  examine  what  the 
singular  line  is  and  what  it  represents. 

The  singular  line  exists  as  a  unique  situation  that 
occurs  while  mapping  the  s-plane  into  the  parameter  plane. 
The  necessary  condition  that  saust  exist  is  that  the 
determinants  of  the  matrices  which  form  equation  (2-20)  must 
be  singular.  tfhen  this  occurs,  a  discontinuity  appears  in 
the  mapping  from  the  s-plane  to  the  parameter  plane.  The 
result  of  this  discontiniut y  is  the  singular  line. 

The  discontinuity  occurs  for  the  singular  value  of  zeta 
and  the  singular  value  of  omega,  and  it  maps  into  the 
parameter  plane  as  a  straight  line  from  plus  infinity  to 
minus  infinity.   The  equation  of  this  line  is 


1 


3 


1 


6=  -  - c 

*    c    c 
1    1 


D     B 
2     2 


C     C 
2     2 


oL 


(3-1) 


The   slope   of   the  line  is  -D  /C  and  the  beta  intercept  is 

11 

-D  /C  . 
1   1 


Until  this  writing  it  was  thought  that  is  a  /alue  of 
alpha  and  beta  from  any  point  along  the  singular  line  was 
placed  in  the  characteristic  equation,  the  system  would  have 
a  pair  of  complex  roots  corresponding  the  the  value  of  the 
constant  zeta-constant  omega  of  the  singular  line.  however 
it  is  shown  in  the  following  examples  that  this  result   does 
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not  always  occur, 


EXAMPLE  III 


Consider  the   following   aircraft   pitch  stabilization 
system  [  6  ]. 


Brk 


^ 

G,[S) 

GJ9 

V 

:J 

\y^\p) 

\ 

a 

Hfl 

Figure  3.1 
AIRCRAFT  PITCH  STABILIZATION  SYSTEM 


where:  Q     is  the  pitch  rate  command 
i 

q     is  the  pitch  rate 

q     is  the  pitch  angle 
o 

*?  is  the  elevator  anale 


and:    G  (s)  is  a  pcle-zero  compensator 

G  (s)  is  a  hydraulic  servo 
2 

G  (3)  is  the  airframe 
H (s)  is  a  rate  gyro 
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The    transfer    functions   of    the   system    components   are 


S-K3C 
S+£ 


K    3180 
1 
G2(s)     = 


2 
s      +    54. 85    +    3201.01 


K    (19.41)  (s    +    3.22) 
2 


Vs)  ■  •    2 

s(s      +    7.58s    +    124.614) 
H(s)     =    s 


The    characteristic    equation    for    the  uncompensated   system   is: 

5  4  3  2 

s       +    62.38s       +    3741.01s       +     (31092.5  +    61732. 6K)s 

+     (398891  +    198779K)s    =    0  (3-2) 

where         X   =    K   K 
1    2 


This   system   is  unstable  for  K>3.27.   It  is  desired  to 

compensate  for  K  =  5  and  to  avoid  the  natural  frequencies  of 

the  airframe,  £0  =  3.24,  and  the  hydraulic  servo,  £0  =  7.049 

n  n 

The  characteristic  equation  of  the  compensated  system  is: 


6  5  -  4 

s   +  (62.38  +£) s   +  (3741.01+  +  62.38@)s 

3 
♦  (31092.5  +  61732. 8K  +  3741.01fl)s   +  (398891  +  198751K<£. 

2 
+  31092.  5A)s   +  (198751FA.+  3988910)s  =  0  (3-3) 


This  characteristic  equation  with  K  =  5  was  used  with 
the  program  PARAM  A  to  obtain  the  parameter  plane  carves, 
figure  3.2,  and  with  the  program  SINGULAR  LINE  to   determine 
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the  existauce  of  singular  line.  Two  singular  lines  were 
found  for  each  zeta  value  of  0.1,  0.2,  0.3,  0.4,  and  0.5, 
and  one  singular  line  was  found  for  each  zeta  value  of  0.6, 
0.7,  0.8,  and  0. 9. 

Zeta  =  0.6  was  chosen  for  compensation  as  it  has  an 
omega  of  12.74  radians  and  would  result  in  a  fairly  fast 
response  without  causing  resonance  at  the  natural 
frequencies  of  the  airframe  and  the  servo.  The  singular 
line  for  zeta  =  0.6  is  the  broken  line  on  the  parameter 
plane  of  figure  3.2. 

Three  pcints  were  selected  on  the  singular  line.  The 
points  were;  A{14.56,50),  3(21.13,150)8  and  C  (27 .69  ,250)  . 
These  values  of  alpha  and  beta  -  where  put  intc  the 
characteristic  -equation  and  the  roots  were  solved  for  on  the 
IBM  360  with  a  polynomial  rootfinder  routine. 

There  were  two  real  roots  and  two  pairs  of  complex 
roots  tor  each  of  the  points.  \ Contrary  to  the  supposition 
that  all  points  along  the  singular  line  will  produce  a  pair 
of  complex  roots  corresponding  to  the  value  of  the 
zeta-omega  of  that  line,  in  this  case  such  roots  were  not 
found.  Ihe  singular  line  with  zeta  =  0.6  and  omega  =  12. 74 
should  have  according  to  prevoius  concepts  result  in  a 
complex  rcot  pair  located  at 

s  =  -7.645  ±  j10. 194  (3-4) 

This  was  not  found  to  be  the  case.  The  complex  roots 
that  were  found  did  not  agree  with  the  singular  line  roots, 
but  they  did  agree  with  the  parameter  plane  roots.  These 
results  are  as  follows: 
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"-scale:  20  units/inch 
/8-scale:  50  units/inch 


Fiqure  3.2 


AIRCRAFT  PITCH  STABILIZATION  SYSTEM  PARAMETER  PLANE 
WITH  ZETA=0.6-OMEGA=12.74  SINGULAR  LINE 
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RQIll 

COMPLEX  ROOTS 

ZETA 

OMEGA 

A 

-8.083±j3.718 
-2.512±J69.00 

.908 
.036 

8.897 
69.04 

B 

-7.944±j6. 987 
-16.66±j64.41 

.751 
.251 

10.58 
66.53 

C 

-7.871±i8.012 
-20.75±j59.36 

.701 
.330 

11.23 

62.87 

Identical  results  were  found  when  values  were  taken 
from  other  singular  lines  in  the  same  relative  area  of  the 
parameter  plane.  All  the  roots  were  found  to  agree  with  the 
parameter  plane,  but  not  with  the  singular  lines. 

The  initial  response  to  these  results  was  that  the 
singular  lines  that  were  found  were  incorrect.  However, 
after  carefully  checking  the  program  siat  were  found  were 
incorrect.  However,  after  carefully  checking  the  program 
SINGULAR  LINE  and  hand  calculating  two  of  the  singular 
lines,  the  program  was  deemed  to  be  correct. 

In  an  attempt  to  find  a  solution  to  this  dilemma  the 
singular  line  examples  given  by  Bowie  were  examined.  In  all 
cases  che  singular  lines  were  located  outside  the  region  of 
the  normal  complex  root  parameter  plane  curves. 

In  the  belief  that  the  constant  zeta  curves  in  figure 
3.2  formed  a  loop  as  omega  increased,  values  of  alpha  and 
beta  believed  to  ce  outside  the  region  were  tried.  The 
straight  lines  connecting  the  zeta  curves  in  the  figures 
shown  in  this  thesis  are  the  result  of  the  DHAW  routine  and 
are  not  part  of  the  parameter  plane.  They  result  when 
intermediate  points  are  not  in  the  area  of  the  plot.  Again 
three  points  on  the  same  singular  line  were  selected.  These 
points  were  D  (76. 93 , 1 000)  ,  E ( 142.57,2000)  and 
F  (208. 22,3000)  .  These  points  are  not  located  on  figure  3.2, 
but   the   location   was   computed   using  the  slope-intercept 
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method.   These  results  are  as  follow: 


POINT     COMPLEX  SOOTS 

D       -7.725±j9.542 
-23.31±;j51.08 

I       -7.688±j9.856 
-23.47ij49.55 

F       -7.675±j9.966 
-23.50±J49.04 


ZETA 


OMEGA 


.629 
.415 

12.23 
56.15 

.615 
.428 

12.50 
54.59 

.610 
.432 

12.58 
54.33 

It  is  evident  from  the  above  that  one  pair  of  complex 
roots,  indicated  by  the  astrisk,  is  approaching  the  value  of 
those  which  shculd  result  due  to  the  singular  line  and  that 
the  value  of  zeta  and  omega  for  that  pair  of  roots  is 
relatively  insensitive  to  large  changes  in  alpha  and  beta. 


In  light  of  these  results  an  examination  of  the 
parameter  plane  in  the  region  of  these  points  was  conducted. 
For  reasons  of  clarity  only  three  values  of  zeta  were 
plotted.  These  were  0.0,  0.6,  and  1.0.  The  resulting 
parameter  plane  is  seen  in  figure  3.3.  As  can  be  seen  the 
constant  zeta  lines  do  not  loop  around  as  prevoiusly 
thought,  but  are  discontinuous  and  the  singular  line 
approaches  asymptotically. 

What  has  happened  is  that  as  omega  was  increased  on  the 
constant  zeta  curve  and  approched  the  singular  point  the 
matrix  which  forms  the  denominator  of  equation  (2-20) 
becomes  very  much  smaller  than  the  numerator  causing  the 
constant  zeta  line  to  go  to  plus  infinity  at  the  singular 
point  asymptotic  to  the  singular  line.  Once  the  value  of 
omega  is  past  the  singular  point  it  reenters  the  parameter 
plane  at  minus  infinity,  asymptotic  to  the  singular  line. 
Since  there  is  only  one  singular  point  for  this  value  of 
zeta  the  constant  zeta  line  will  continue  through  the 
parameter   plane   approaching   infinity   once   more  as  omega 
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approaches  infinity. 

The  reason  that  this  occurs  is  evident  when  the  mapping 

of  the  constant  zeta  line  form  the  s-plane  to  the   parameter 

plane  is  observed.   The  constant  zeta  line  in  the  s-plane  is 

continuous  fcria  CO  =0  to  CO  =00,  therefore  the  mapping  of   this 

n       n 

line   into   the   parameter   plane   must   also  be  continuous. 

However  there  is  a  discontinuity  at  the  singular  point  where 

the   constant   zeta   curve   in   the   paramter   plane  goes  to 

infinity  for  a  finite  value  of   omega.    The   singular   line 

provides   a  continuous  mapping  of  the  constant  zeta  curve  as 

it  traverses  from  plus  infinity   to   minus   infinity.    This 

line   exists   at  infinity,  but  not  in  the  parameter  plane  as 

it  is  known. 

In  an  effort  to  examine  the  change  in  alpha  and  beta  as 
omega  was  increased  a  printout  of  PA SAM  A  was  made.  However 
due  to  the  iteration  of  omega  in  PASAM  A  it  was  not  possible 
to  obtain  sufficient  data  around  the  singular  point.  In 
PARAM  A  for  an  omega  of  about  12  radians  the  iterative  step 
is  about  0.9  radians.  Therefore  the  section  of  PARAM  A 
which  computes  the  constant  zeta  curves  was  modified  to  vary 
omega  over  1  radian  in  steps  of  0.001  radians  starting  0.5 
radians  below  the  singular  point  and  increasing  to  0.5 
radians  above  the  singular  point. 

In  order  to  see  how  small  increases  in  omega  affect  the 
values  of  alpha  and  beta  some  of  the  points  are  shown: 


OMEGA 

12.242 
12.342 
12.442 
12.542 
12.642 
12.741 
12.743 
12.842 
12.942 
13.042 


ALPHA 


BETA 


88.56 

1297.87 

111.31 

1644.01 

149.21 

2220.85 

226.06 

339  1  .11 

451.58 

6826. 15 

144447. 

2200360. 

66473.4 

-1012672. 

-454. 96 

-6980.32 

-230.44 

-3564.49 

-154.32 

-2405.44 
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Figure  3.3 


PARAMETER  PLANE  FOR  AIRCRAFT  PITCH  STABILIZATION  SYSTEM 

ZETA  =0.6 
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As  can  be  seen  for  large  changes  in  alpha  and  beta  in  the 
area  of  the  parameter  plane  where  the  singular  line  is 
asymptotic  to  the  constant  zeta  curves  results  in  small 
changes  of  omega.  Since  oaega  is  very  insensitive  to 
changes  on  alpha  and  beta  in  this  region  the  concept  of  the 
singular  line  is  valid  from  the  standpoint  of  the  engineer. 
It  is  possible  to  vary  a  parameter  along  the  singular  line 
in  this  region  and  expect  to  maintain  a  constant  value  of 
zeta  and  omega. 

EXAMPLS  IV 

In  light  of  the  unexpected  results  obtained  in  example 
III  it  was  decided  to  carefully  reexamine  example  II  that 
was  given  in  the  previous  chapter.  The  parameter  plane  was 
expanded  around  the  area  of  the  origin  where  the  singular 
lines  pass  through  the  complex  root  portion  of  the  parameter 
plane. 

Figure  3.4  shows  this  area  of  the  parameter  plane.  The 
constant  7eta  curves  are  for  values  of  0.0,  0.2,  0.5,  0.6, 
0.8,  and  1.0.  There  are  two  singular  lines  for  each  zeta 
value  of  0.2,  0.5,  0.6,  and  0.8  which  are  shown  as  dashed 
lines.  There  are  two  discontinuities  in  each  zeta  curve. 
Although  it  is  not  immediately  obvious  due  to  the  limited 
area  of  the  plot,  all  the  constant  zeta  curves  originate  in 
the  second  guadrant,  and  become  infinite  at  the  first 
dicontiuity  and  then  reappear  at  infinity  in  the  fourth 
guadrant.  The  curves,  for  values  of  zeta  up  to  and 
including  0.5  then  follow  a  path  just  to  the  left  of  the 
origin  and  enter  the  first  quadrant  where  they  again  become 
discontinuous  and  infinite.  The  zeta  curves  for  values  of 
zeta  greater  than  0.5  come  from  infinity  in  the  fourth 
guardant  and  go  to  the  left  of  the  origin  where  they  turn 
into   the    third    quadrant    where    they    again    become 
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discontinuous  and  infinite.  The  curves  for  zeta  values  of 
0.5  and  less  reappear  at  infinity  in  the  third  guadrant  then 
go  to  infinity  as  a  continuous  function  in  the  first 
guadrant.  The  curves  for  zeta  values  greater  than  0.5 
reappear  at  infinity  in  the  first  quadrant,  sweep  toward  trie 
origin  then  tack  through  the  first  quadrant  to  infinity 

The  region  of  stability  in  the  parameter  plane  is 
bounded  by  the  zeta=0.0  curve,  therefore  the  singular  lines 
of  interest  lie  in  the  first  guadrant.  As  can  be  seen  the 
zeta=0.5  and  zeta=0.6  singular  lines  nearly  coincide  with 
the  corresponding  zeta  curves  in  the  parameter  plane  and  tne 
zeta=0.2  and  zeta=0.8  lines  are  being  asymptotically 
approached  by  the  corresponding  parameter  plane  zeta  curves. 

As  in  example  III  the  modified  section  of  PARAK  A  »as 
used  to  abtain  a  printout  of  alpha  and  beta  for  a  1  radian 
change  in  omega  as  the  constant  zeta  curve  passed  through 
the  singular  line  value.  In  this  example  the  constant  zeta 
curves  in  the  first  guadrant  approached  their  singular  value 
very  rapidly. 

For  purposes  of  illustration  points  have  been  placed  on 
the  constant  zeta  curves  in  figure  3.4  at  the  value  at  which 
omega  is  within  0.5  radians  of  the  singular  point  which  is 
at  infinity.   These  points  are: 


POINT 


ZETA 


ALPHA 


BETA 


A 

0.20 

4.078 

1.30 

E 

0.50 

-0.542 

0.59 

C 

0.60 

0.571 

1  .60 

c 

0.80 

4.754 

6.62 
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Figure  3.4 


SINGULAR  LINE  THROUGH  THE  PARAMETER  PLANE 
FOURTH  ORDER  POLYNOMIAL 
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As  can  be  readily  seen  these  points  are  within  5 
percent  cf  the  singular  value  of  omega  and  if  one  continues 
along  the  zeta  curves  it  will  be  found,  that  for  practrical 
design  purposes,  the  changes  in  omega  are  inconsequential. 


It  has  been  shown  that  a  true  singular  line  in  the 
original  concept  aces  not  exist.  In  examples  III  and  IV  it 
is  shown  that,  if  a  designer  carefully  examines  the  singular 
lines  and  their  relationship  to  the  constant  zeta  curves  of 
the  parameter  plane  the  singular  concept  can  be  used  for 
practical  design  applications. 
EXAHPLE  V 

This  example  is  the  original  system  in  which  Bowie 
found  that  singular  lines  existed  in  the  parameter  plane. 
The  system  which  was  used  has  a  characteristic  equation 
which  is  the  following  sixth  order  polynomial: 

6       5  4         3  2 

s   +  80s   +  (2001  +  1600)  s   +  84Gv<s   -f-  (1 600cU400£)  s 

+  1600£s  *  1600$  =  0      (3-5) 

The  parameter  plane  with  constant  zeta  curves  with  the 
values  of  0.0,  0.2,  0.4,  0.6,  0.8,  and  0.9  is  seen  in 
figures  3.5  and  3.6.  Figure  3.6  is  the  area  around  the 
origin  in  figure  3.5.  There  were  two  singular  lines  tor 
each  of  the  zeta  values  0.2,  0.4,  0.6,  0.8,  and  0.9  and 
these  are  shown  on  figures  3.5  and  3.6  as  dashed  lines. 

As  can  be  seen  the  parameter  plane  has  no  visible 
discontinuities  as  omega  increases,  however  due  to  the  fact 
that  two  singular  lines  exist  for  each  value  of  zeta  there 
roust  be  two  discontinuities  in  each  zeta  curve.  Examination 
of  the  singular  lines  and  the  constant  zeta   curves   of   the 
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parameter  plane  indicates  chat  the  singular  lines  are 
tangent  tc  the  corresponding  zeta  curve  at  one  point.  When 
the  values  of  zeta  and  omega  for  these  points  were  examined 
it  was  found  that  the  tangent  points  occur  at  the  singular 
points  in  the  constant  zeta  curves. 

Unlike  examples  III  and  IV  where  the  mapping  from  the 
s-plane  approached  infinity  as  the  singular  pcint  was 
approached  this  example  maps  as  though  no  discontinuity 
existed.  When  the  discontinuity  is  reached  the  constant 
zeta  curve  immediately  goes  to  plus  infinity  tangent  to  the 
curve  at  the  point  of  discontinuity  and  returns  from  minus 
infinity  tangent  to  the  constant  zeta  curve  and  continues 
mapping  to  infinity  as  omega  approaches  infinity. 

This  phemonenon  does  not  reveal  itself  in  a  normal 
parameter  plane  study  for  two  reasons.  The  first  of  these 
is  the  previousJ.y  mentioned  size  of  iteration  in  the  PARAM  A 
program  and  the  second  reason  is  that  a  normal  practice  in 
computer  programing  is  to  limit  the  numerator  of  a  term  to  a 
minimum  size  to  prevent  an  overflow  or  underflow  in  the 
computer  when  a  number  beyond  the  computer's  capablity  is 
reached.  This  practice  was  followed  in  the  writing  of  PAKAM 
A  so  if  the  singular  point  was  reached  on  an  iterative  step 
the  program  would  have  recognized  this  as  a  forbidden  number 
and  continued  to  the  next  value  of  omega,  thereby  bypassing 
the  singular  point  and  continuing  with  no  indication  that 
this  condition  existed.  If  Bowie  had  not  oeen  interested  in 
the  area  outside  the  normal  complex  root  area  of  this 
example  he  would  have  had  no  indication  that  singular  lines 
existed. 
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«-scale:  50  units/inch 
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Figure  3.5 


SINGULAR  LINE  IN  THE  PARAMETER  PLANE 
EXAMPLE  V 
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"-scale:  10  units/inch 
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Figure  3.6 


SINGULA!;  LIMES  IN  THE  PARAMETER  PLANE 
EXAMPLE  V  -  EXPANDED 
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In  this  example  the  singular  lines  are  true  singular 
lines  as  originally  conceived.  If  a  value  of  alpha  and  beta 
is  selected  at  any  point  along  these  lines  a  pair  of  complex 
roots  corresponding  to  the  constant  zeta-constant  omega  of 
that  line  will  result.  To  show  this  fact  six  points  have 
been  selected  in  the  region  of  singular  line  on  figure  3.5. 
As  can  be  seen  by  the  following  the  expected  complex  roots 
were  obtained. 


'IliJ 

HIE 

A 

0.80 
0.40 

B 

0.90 
0.60 

C 

0.80 
0.60 

D 

0.90 
0.60 

E 

0.60 
0.40 

F 

0.60 
0.80 

SINGULAR  POINT 


COMPLEX  ROOTS 
EXPECTED  OBTAINED 


23 

47 

20 
31 


55 
41 

88 
53 


3.396 
31.53 

3.831 
31.53 

2.53  7 
47.41 

2.537 
23.55 


•16.44+ 
18.96± 

■18.79± 
13.92± 

2.717± 

18.92± 

•3.450± 
■18.92± 

1.522± 

18.96± 

■1.522± 

•18.44± 


114.  13 
J43.15 

19.101 
J25.22 

12.078 
;J25.22 

10.728 
J25.22 

12.030 
"J43.45 

12.030 
J14.13 


18 
18 

18 

18 


42±1 
95±j 

7  9±1 
9  1±j 

2.721±j 
18.91±j 

3.465±1 
18.90±j 

•1.524±j 
•18.94±3 

1.521+j 

18.47±j 


14.  15 

43.4  4 

9.  121 

25.23 

2.081 
25.22 

0.729 
25.23 

^.032 
43.4  3 

2.029 
14.15 


In  summary,  the  singular  lines  result  when  a  constant 
zeta  or  constant  omega  curve  is  mapped  from  the  s-plane  to 
the  parameter  plane  and  a  discontinuity  occurs.  At  this 
discontinuity  the  parameter  plane  mapping  goes  to  infinity, 
either  plus  or  minus,  and  returns  from  the  other  side  of 
infinity.  The  singular  line  simply  provides  a  mapping 
through  this  region  to  maintain  a  continuous  mapping  from 
the  s-plane  to  the  parameter  plane.  In  cases,  such  as 
examples  III  and  IV,  where  the  mapping  first  approaches 
infintiy  as  the  singular  point  is  approached  the  singular 
line  concept  is  a  valid  engineering  aid  over  a  limited  range 
of  the  singular  line.  However  in  cases  such  as  example  V 
the  singular  line  may  be  used  over  its  entire   range   as   an 
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excellent    design   aid. 
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IV.   THE  SINGULAR  SURFACE  IN  THE  PARAMETER  SPACE 

A.   PARAMETER  SPACE 

The  parameter  space  is  an  extension  of  the  parameter 
plane  when  there  are  three  or  more  variable  parameters  [7], 
[8],  [9].  Once  the  concept  of  the  parameter  plane  is 
understood  it  is  relatively  easy  to  extend  the  concept  to 
three  or  more  variable  parameters. 

The  easiest  method  to  generate  a  three  dimensional 
parameter  space  is  to  set  the  third  variable  parmeter  the 
third  variable  parameter  to  a  known  value  and  then  generate 
the  two  dimensional  parameter  plane.  The  process  is  then 
repeated  for  additional  values  of  the  third  parameter.  The 
result  is  a  series  .of  parameter  planes  which  will  form  a 
parameter  space. 

EXAMPLE  VI 

Consider  the  following  system: 


■>> 

S+  a 

s+£ 

s+20 
s+30 

K 

v 

J 

s  (s+1)   (s+2) 

\ 

Figure  4. 1 
THIRD  ORDER  SYSTEM  WITH  DUAL  CASCADE  COMPENSATION 
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The    characteristic    equation    for    this    system    i.«.: 

5  4  3  2 

s      +     (33    +^)s      +     (9  2   +    32£)  s      +     (60    +    92£  +    K)  s 

+     (2  OK    +<*K    +    60)s    +    Kot=0  (4-1) 

Figure  4.2  shows  the  parameter  space  for  the  constant 
zeta=0.5  and  constant  omega=1.7  and  the  values  of  K  from  10 
to  15.  It  was  constructed  by  fixing  the  value  of  K  and 
computing  the  parameter  plane  using  PARAM  A.  The  curves  were 
then  plotted  onto  a  single  plot  to  show  the  three 
dimensions. 

Both  the  constant  zeta  curves  and  the  constant  omega 
curves  form  a  surface  and  the  interconnection  of  these 
surfaces  forms  a  curved  line  through  the  parameter  space. 
Any  point  on  this  line  will  generate  a  pair  of  complex  roots 
corresponding  to  the  s-plane  point  of  zeta=0.5  and 
omega=1.7.  For  these  values  of  zeta  and  omega  the  complex 
roots  are: 


-0.850  ±  J1.470 


(4-2) 


By  selecting  three  points  on  this  line  this  is  shewn  to 
be  true.  The  points  ere  A  (0.16,  2.6,  10.0).  B(0.347,  3.0r 
12.0),  and  C  (0.52,  3.63,  15.0).  These  points  were  taken 
from  the  parameter  plane  using  an  engineer's  scale. 


POINT 

COMPLEX  ROOTS 

SETA 

Of*  EG  A 

A 

-0.8t7±j1.478 

.497 

1.703 

B 

-0.8i±9±j  1  .  475 

.498 

1.701 

C 

-0.851±j1.469 

.501 

1.698 
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Figure  4.2 
PARAMETER  SPACE  WITH  ZETA=0.5  AND  0MEGA=1.7  SURFACES 
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B.   SINGULAR  SURFACES 

The  parameter  space  is  continuous  throughout  the  range 
of  the  third  variable  parameter.  If  for  any  value  of  the 
third  parameter  a  discontinuity  appears  in  the  parameter 
plane  it  will  also  appear  in  the  parameter  space  as  a 
discontinuity.  As  stated  in  chapter  III  this  discontinuity 
will  result  in  a  singular  line.  As  the  third  parameter  is 
varied  the  discontinuity  will  result  in  sdngular  lines  for 
the  same  value  of  zeta  and  omega.  The  slope  and  the  axis 
intercept  points  will  change  slighty  so  the  result  is  a 
smooth  surface  representing  a  warped  hypersurface  traversing 
through  the  parameter  space. 

Figure  4.3  is  a  photograph  of  a  singular  surface  model. 
The  singular  surface  is  from  the  system  in  example  VI.  the 
singular  values  are  zeta=0.5  and  omega=12.04  and  the  third 
variable  K  is  varied  from  5.0  to  20.0.  As  can  be  seen  the 
surface  is  a  smooth  warped  hyperplane.  The  model  was 
constructed  by  obtaining  the  singular  lines  for  the  same 
singular  value  at  each  value  of  K,  transferring  them  to 
cardboard,  cutting  them  and  accentting  them  with  a 
contrasting  color.  The  photography  was  done  by  the  Photo 
Division,  Education  Media  Department,  Naval  Postgraduate 
School. 

Figure  4.4  is  a  drawing  of  the  same  singular  surface  as 
figure  4.3.  The  difference  between  the  drawing  and  the 
photograph  is  that  in  the  photograph  the  K  axis  starts  at 
K=5  in  the  foreground  continues  back  through  to  K=20,  while 
in  the  drawing  K=5  is  the  most  distant  singular  line  and  the 
singular  lines  progress  forward  until  the  value  of  K=20  is 
reached. 
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Figure    4. J 
SINGULIS    SURFACi    MODEL 
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Figure  4.4 
SINGULAR  SURFACE  DRAWING 
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The  relationship  of  the  singular  surface  to  the 
parameter  space  is  exactly  the  same  as  that  of  the  singular 
line  to  the  parameter  plane.  As  with  the  singular  line,  the 
singular  surface  provides  continuity  at  infinity.  The 
designer  must  Le  aware  of  the  type  of  discontinuity  which  is 
producing  the  surface.  If  the  discontinuities  occur  like 
the  ones  in  example  III  and  IV  then  the  principles  -of 
singular  lines  may  only  be  applied  in  limited  areas  of  the 
surface.  If  the  singularities  occur  like  those  in  example 
IV  then  any  point  on  the  singular  surface  will  give  a  pair 
of  complex  roots  corresponding  to  the  constant  zeta-constant 
omega  of  the  surface. 
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3Ls  compensation  techniques  using  the  singular  line 

The  major  advantage  of  using  singular  lines  for 
compensation  is  that  as  one  moves  up  and  down  the  singular 
line,  witnin  previously  mentioned  limits,  one  pair  of 
complex  roots  remains  in  the  same  locations  in  the  s-plane. 
While  this  pair  of  roots  is  stationary  other  root  locations 
can  be  varied  as  the  location  on  the  singular  line  is  moved 
through  other  singular  lines,  real  root  lines,  and  other 
constant  zeta  parameter  plane  curves. 

There  many  excellent  texts  available  which  go  into 
great  detail  on  how  root  location  will  affect  system 
response,  overshoot,  damping  and  so  forth  [1],  [10].  Once 
the  basic  principles  on  how  root  locations  affect  a  system 
are  understood  it  is  possible  to  use  singular  line  concepts 
to  obtain  proper  root  locations. 
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EXAMPLE  VIII 

Consider    the   pitch   stabilization   system   given   in 
example  III. 


a 


Hfl 


Figure  5, 1 
AIRCRAFT    PITCH  STABILIZATION  SYSTEM 

The  characteristic  equation  for  the  system  is; 


6  5  4 

s      +     (62.33    +^)  s       +     (3741.01+    +    62.33p)s 

3 
+     (31092.5    +    61732. 8K    +    374  1.010)  s      +     (398891    +    193751ft* 


+    31092.  5£)  s      +     (198751KOC+    398891£)s    =    0 


(5-1) 


It  was  shofen  in  example  III  that  singular  lines  exist  for 
this  system,  but  they  are  the  result  of  discontinuities  at 
infinity.  Therefore  in  order  to  effectively  use  the 
singular  line  concept  it  is  necessary  to  examine  the 
singular  lines  to  determine  the  areas  in  which  singular 
results  could  be  expected. 

In   the  ca.se   of  this  problem  it  was  decided  to  use  the 
singular  line  with  zeta=0.6  aid  oraega=12.74   as   a   dominant 
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root  pair  to  give  fast  system  response  and  avoid  the  natural 

freguenci.es  of  the  servo  and  the  airframe.  Figure  5.2  shows 

the   singular   lines   tor   this  system  for  values  of  zeta  of 
0.4,  0.5,  0.6,  0.7,  0.8,  and  0.9. 

By  moving  up  the  zeta=0.6,  omega=12.74  singular  line 
above  the  region  where  beta=3000  the  value  of  zeta  and  omega 
should  remain  relatively  constant.  Two  points  have  been 
selected  to  attempt  compensation,  A (273. 87,  4000)  and 
B(306.69,  4500).  These  values  have  been  put  into  the 
characteristic  equation  and  the  system  was  simulated  using  a 
step  input  of  0.2  radians  per  second. 

Figure  5.3  is  the  response  for  the  point  A  values  of 
alpha  and  beta  and  figure  5.4  is  the  response  for  the  point 
B  values  of  alpha  and  beta.  As  can  be  seen  the  system 
response  is  essentially-  the  same  for  both  points.  The 
response  from  the  point  A  values  has  a  peak  overshoot  of 
1.08  and  the  response  from  the  point  3  values  has  a  peak 
overshoot  of  1.1.  The  steady  state  output  is  1.66  degrees 
in  both  cases. 

This  demonstrates  one  of  the  main  advantages  of  the 
singular  line  concept.  As  two  parameters  are  varJLed  the 
system  response  remains  insensitive  to  the  change. 

If  this  parameter  variation  were  to  take  place  inside 
the  normal  ccmplex  root  region  of  the  parameter  plane  the 
shift  of  root  locations  would  result  in  a  change  of  output 
response.  Figure  5.5  is  the  parameter  plane  of  the  normal 
complex  root  area  of  the  airdraft  pitch  stabilization 
system.  Between  point  A  and  point  C  the  linear  change  in 
beta  is  200  and  the  change  in  alpha  is  13.13. 

There  are  two  pairs  of  complex  roots  and  two  real  roots 
in   this   section   of   the  parameter  plane.   One  of  the  real 
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roots  is  cancelled  at  the  origin  of  the  s-plane  by  the  zero 
produced  by  the  rate  gyro  in  the  feedback  loop.  The  other 
real  root  location  is  past  s=-75  in  these  cases  and  will 
only  have  a  slight  affect  on  the  steady  state  value  of  the 
output.  The  complex  root  locations  and  their  corresponding 
values  of  zeta  and  omega  are: 

POINT    co^jl^Z.   £OOTS    zeta    2M£A 

A       -8.0d3± 13-713     .908     8.897 
-2.512±j69.00     .036     69.04 

B       -7.944±j6.987     .751     10.58 
-16. 66±:J64.41     .251     66.53 

C       -7.871+18.012     .701     11.23 
-20.75±J59.36     .330     62.87 


At  point  A  the  dominant  root  is  located  at  s  =  -2.512  ± 
J69.00.  At  point  C  the  dominant  root  is  located  at  s  = 
-7.871  ±  J8.012.  It  can  be  seen  that  at  the  point  A  values 
the  system  response  will  be  almost  undamped  with  a  peak 
overshoot  of  about  1.9  and  a  natural  frequency  of  about  69 
radians,  while  at  the  point  C  values  the  system  will  be  well 
damped  with  a  natural  frequency  of  about  11  radians. 
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"-scale:  200  units/inch 
£-scaie:  1000  units/inch 


Figure  5.2 


SINGULAR  LINES  FOR 
AIRCRAFT  PITCH  STABILIZATION  SYSTEM 
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Figure  5.3 
STEP  RESPONSE  FOR  POINT  A 
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Figure  5. '4 
STEP  RESPONSE  FOR  POINT  B 
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a-scale:  20  units/inch 
/?-scale:  50  units/inch 

Figure  5.5 
AIRCRAFT  PITCH  STABILIZATION  SYSTEM  PARAMETER  PLANE 
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VI.   CONCLUSIONS  AND  RECOMMENDATIONS  FOR  FURTHER  SIUDY 


It  has  been  shown  that  some  of  the  original  concepts  of 
the  singular  line  were  in  error  due  to  previously  unknown 
properties  ox  the  singular  line  and  its  relationship  to  the 
parameter  plane.  The  main  reason  these  properties  were  not 
found  earlier  was  due  tmainly  to  standard  computer 
programming  practices  and  limitations. 


Although  it  has  been  found  that  some  limitations  exist 
when  using  singular  line  concepts  it  has  been  shown  that  the 
basic  concepts  of  singular  lines  can  be  used  by  the  design 
engineer  to  easily  obtain  desired  system  response  when  more 
than  one  variable  parameter  is  present  in  the  system. 

An  example  has  been  presented  which  shows  that  the 
singular  line  concept  is  an  excellent  method  to  obtain 
desired  system  response.  The  method  will  allow  the  designer 
to  select  either  dominant  root  locations  or,  if  it  is 
desired,  nondominant  root  locations.  the  root  locations 
when  selected  from  the  singular  line  region  of  the  parameter 
plane  will  remain  fixed  in  their  location  as  long  as  the 
alpha  and  beta  parameters  are  restricted  to  the  singular 
line. 


For  further  study  it  is  highly  recommended  that  the 
case  where  true  singular  lines  exist  such  as  in  example  V  be 
thoughly  studied  to  obtain  a  full  understanding  of  why  the 
singular  point  does  not  occur  at  infinity  in  the  parameter 
plane.  Once  this  is  accomplished  devise  a  method  to  force 
this  type  of  condition  at  a  desired  zeta-omega  value.  If 
this   can   be   accomplished   it   would  allow  the  designer  to 
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select  a  root  location  in  the  s-plane  and  use  singular   line 
concepts  for  easy  design  and  compensation. 
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APPENDIX  A 


PARAM  A 


This  is  a  modified  version  of  the  PARAM  A  program 
written  by  Nutting  [5  J.  Originally  written  in  FORTRAN  for 
the  CDC  1604  Computer,  it  has  been  modified  to  FORTRAN  IV 
for  the  IBM  360/67  at  Naval  Postgraduate  School.  It  has 
also  been  modified  to  allow  the  user  to  select  the  size  of 
the  output  plot  and  to  plot  singular  lines  in  the  parameter 
plane. 

However   PARAM   A   will   not  determine  the  existance  of 

singular  lines.   It  is  necessary  for  the  user   to   know   the 

value   of  c  and  w   of  a  singular  line  before  this  opticn  can 

n 

be  used.   These  values  can  best  be  determined  by  use  of   the 

program   SINGULAR   LINE.   The  user  then  can  either  put  these 

values  into  PARAM  A  or  use  the  values  to  draw   the   singular 

line  on  the  PAEAM  A  output  plot. 
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APPENDIX  B 


SINGULAR  LINE 


The  purpose  of  the  program  SINGULAfi  LINE  is  to 
determine  the  existance  of  singular  lines  in  the  parameter 
plane.  Ihe  program  is  written  in  FORTRAN  IV  for  the  IBM 
360/67  computer  and  the  CALCOMP  plotter  at  Naval  Postgrauate 
School. 

The  values  read  into  the  program  are  the  coefficients 
of  the  characteristic  equation,  the  values  of  zeta  to  be 
examined,  and  plotting  information.  The  program  contains 
three  subroutines  which  are  located  in  the  source  library  of 
the  IBM  360/67.  These  are  POLRT  which  has  the  ability  to 
solve  up  to  a  36  order  polynomial,  DRAW  which  outputs  the 
plot  to  the  CALCOMP  plotter,  and  BCNV  which  converts  numeric 
symbols  to  alpha-numeric  for  use  in  DRAW. 

The    program   takes   the   values   of   zeta   and   first 

generates  the  Chebyshev  function  of   the   second   kind.    It 

then   generates   the  coefficients  for  the  omega  terms  in  tne 

polynomial  expansion  of  B  C   -  3  C   =0.   Ihis  generation  is 

12     2  1 

accomplished  by  using  the  coefficients  of  the  alpha  and  beta 

terms  in  the  characteristic  equation,  the  value  of  zeta   and 

the   Chebyshev   function.    The  resulting  omega  coefficients 

are  now  input  into  the  subroutine  POLRT. 

POLRT  solves  for  the  roots  of  the  omega  polynomial  and 
outputs  the  results  to  the  program.  These  roots  are  now 
examined  for  positive  real  values.  If  positive  real  roots 
are   present   there   is   a   possibility   of   singular   lines 
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existing  for  that  value  of  zeta.  If  there  are  not  positive 
real  roots,  there  are  no  singular  lines.  This  information 
is  printed  out  and  the  next  value  of  zeta  is  examined. 

For  the  positive  real  roots  B,B,C,C,D,  and   D 

12    12    1         2 

of   equation   (2-27)   are   formed   and  used  to  insure  -chat  a 

singular  line  does  exist.   If  a   singular   line   does   exist 

then   B  ,   C  ,   and  D   are  used  to  generate  tne  points  along 
1     1         1 

the  line  by  using  equation  (2-28).    These   points   are   now 

input   into  the  subroutine  DRAW  and  are  plotted  and  labeled. 

The  printed  output  includes  the  value  of  zeta  and 
omega,  the  alpha  and  beta  intercept  points,  the  slope,  and 
the  label.  If  there  are  less  than  two  points  generated 
within  the  area  of  the  plot  there  will  be  no  plot  of  that 
line.  This  information  is  indicated  on  the  printed  output 
and  the  user  can  either  change  the  plot  scales  or  use  the 
output  line  information  to  construct  a  singular  line  plot. 

This  procedure  is  repeated  for  each  value  of  zeta, 
plotting  all  the  lines  on  a  single  plot.  Multiple  runs  for 
different  characteristic  equations  are  possible  by  using 
additional  data  decks.  Each  equation  will  have  its  own 
output  plot. 
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